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' Abstract 

' , Requiring that the baryon number of the universe be generated by anomalous elec- 
Mh. 

^ ■ troweak interactions places strong constraints on the minimal supersymmetric standard 

^ I model. In particular, the electric dipole moment of the neutron must be greater than 

' -97 

^ ■ 10 e-cm. Improvement of the current experimental bound on the neutron's electric 

' dipole moment by one order of magnitude would constrain the lightest chargino to be 
lighter than 88 GeV, and the the lightest neutralino to be lighter than 44 GeV. In exten- 
sions of this model with gauge singlet superfields all of these bounds are eliminated. 
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Despite the success of the standard model of weak and strong interactions, we stiU 
remain ignorant of the mechanism of electroweak symmetry breaking; we are just beginning 
to probe this sector directly through Higgs searches at LEP. We are also still ignorant of 
the origins of the CP violation observed in the kaon system. Unfortunately, the minimal 
standard model accounts for all solid experimental results observed to date, and we have 
few experimental constraints on the symmetry breaking or CP violating aspects of the 
theory. There is a strong cosmological argument that the minimal standard model cannot 
be the whole particle physics story. The baryon to entropy ratio of the universe is (0.4 — 
1.0)10"^*^ and explaining this observable requires baryon number violation coupled 
with out-of-equilibrium CP violation in the early universe [^. While the standard model 
does have CP violation, the effects of the CP violating phase in the Kobayashi-Maskawa 
matrix are too suppressed by small masses and mixing angles in order to account for the 
observed baryonic asymmetry of the universe (BAU) 0. The standard model also contains 
anomalous baryon number violating interactions [0, which should be rapid enough at 
high temperatures to affect cosmology [0-0. Furthermore, the standard model can satisfy 
the out of equilibrium condition for baryogenesis if the phase transition is first order, 
proceeding via nucleation and expansion of bubbles of the broken phase . Unfortunately 
it is necessary for the vev of the Higgs field after the transition to be large in order to avoid 
washing out any baryon number created during the transition P, p!0| - |r^ . In the minimal 
standard model this requirement cannot be satisfied unless the Higgs is lighter than ~ 35 



GeV [|T2[, which conflicts with current experimental bounds. Thus we should look beyond 
the standard model in order to discover the origin of the BAU. However we may not need 
to look very far, as possible mechanisms for baryogenesis have been suggested in several 



reasonable extensions such as axion models with additional light scalar doublets the 
singlet Majoron model the two Higgs model [pr7| - p!9| , the supersymmetric standard 
model [p!8| , p!9[| , extended supersymmetric models [ pO| and left-right symmetric models ||21 



Requiring that sufficient BAU be generated during the weak transition, and that it not be 
washed out afterwards can give us new information about the CP violating and symmetry 
breaking sectors of the weak interactions, allowing us to rule out some models (such as 



^ Shaposhnikov has suggested two conceivable ways to enhance the CP violation in the standard 
model at high temperature Q ; the first mechanism, dynamical high temperature spontaneous 
CP violation, is contradicted by non-perturbative computation and the second mechanism, 
reflection of baryon number from expanding bubble walls, according to our estimates cannot 
provide a large enough asymmetry H. 
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the minimal standard model), and to constrain parameters in others. It is the aim of 
this letter to use baryogenesis to find new constraints on supersymmetric models. We are 
able to severely constrain almost all unknown parameters of the minimal model, while 
considerable freedom remains for models with additional gauge singlets. 

First let us consider baryogenesis in the minimal supersymmetric standard model 
(MSSM) p2|. Several authors have claimed that this model is ruled out for baryogenesis 
because the bound on the lightest scalar mass in this model is the same as in the minimal 
standard model [^|T3[. However Myint finds that the bounds on the scalar mass 
are relaxed somewhat, to 64 GeV, due to top quark and squark corrections to the high 
temperature effective potential when the top quark is heavy and the squark masses are 
not too heavy. The bound is relaxed because top squarks can play the role envisioned by 
Anderson and Hall for a gauge singlet scalar, whose coupling to the Higgs doublet increases 
the strength of the transition [^. The upper bound on the Higgs mass was computed 
in one loop perturbation theory and higher order corrections to the gauge propagator 
will reduce this number by a factor of about a/2/3 [^, to ~ 50 GeV. Furthermore this 
bound will receive corrections proportional to {mn /fnz)'^ from two loop corrections. Thus 
avoiding baryon number washout is difficult in the MSSM, but not impossible. To increase 
the upper bound on the Higgs mass as much as possible, one should take the top quark to 
be heavy ( > 150 GeV), the squark masses not much heavier than 150 GeV, the trilinear 
soft supersymmetry breaking terms ( "A-terms" ) small, and the parameter tan (3 (the ratio 
of the two Higgs vevs) less than 1.7. Then, in improved one loop perturbation theory, 
avoiding baryon number washout requires the lightest scalar mass to be lighter than ~ 50 
GeV, which is not in conffict with current bounds provided that higher order corrections 
do not decrease the mass bound. 

The MSSM has two possible sources of CP violation which are absent in the min- 



imal standard model [24]; a combination of these will be constrained by the BAU. The 
interactions in this model are given by 

\U\uQH + DXdQH' + EXeLH' + 1^1 e-^'^^ra'l ^ 

_ _ _ (1) 

+ m3/2 [\A\e''^^{UiuQH + DioQE' + EiELH') + \iiB\HH']^ . 

Note that unlike the generic two Higgs doublet model discussed in refs. [jT^|T^, the Higgs 
potential does not contain any CP violating phases at tree level — the phases 4>a and 
occur only in interactions involving the super-partners of the ordinary particles. We find 
that the contribution of the phase (j)A to the BAU is small, and so we focus on the effects 
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of the phase (pB, which appears in the mass matrices of the supersymmetric partners of 
the gauge bosons and the Higgs scalars (the "mos"). 

Our strategy for calculating the BAU in the MSSM is as follows. We will fix the top 
mass at 170 GeV, the lightest scalar mass at 48 GeV, the soft supersymmetry breaking top 
squark masses at 150 GeV, the A-terms at zero, and /3 at 0.85 (tan/3 = 1.14), since these 
parameters are already constrained to be near these values and in any case allowing 
them to vary will only affect the BAU by 0{1). Using these values we then calculate the 
critical temperature, the effective potential at the critical temperature, and the shape of 
the bubble walls. (Note that these phase transition quantities are not very sensitive to 
other supersymmetric parameters such as ino masses.) We use the improved one loop 
approximation for the effective potential, including the order corrections to the gauge 
propagators, and neglecting the contribution of Higgs doublet loops. When the mass of 
the lightest scalar is far below the gauge boson masses this approximation is reliable at the 
critical temperature around the symmetry breaking minimum and in the vicinity of the 
symmetric minimum for scalar field expectation values larger than 0{g^\J^). Although 
perturbation theory is not valid for a calculation of the effective potential between the 
two minima, which will affect the width and shape of the phase boundary, fortunately 
our calculation of the BAU will turn out to be insensitive to the detailed shape of this 
boundary, provided it is much thicker than 0{1)/T. The one loop estimate gives the 
width of the phase boundary to be 11/T at the temperature where the two minima are 
degenerate, so we will assume the phase boundary is thick. The one loop effective potential 
indicates that the ratio of the expectation values of the two Higgs remains constant during 
the transition, so we take 



where {H') and {H) are the expectation values at the critical temperature in the symmetry 
breaking minimum. The temperature of the transition! is 59 GeV, and at this temperature 
the minimum of the effective potential occurs at -ff = 63 GeV and H' = 53 GeV. 

^ Our definition of tlie transition temperature is tlie spinodal point wliere tfie local minimum 
at the origin vanislies. Tliis temperature is slightly lower than the one loop estimates of the 
temperature at which the transition actually occurs p^ , but these estimates require knowing the 
effective potential in a region where it is not calculable perturbatively. We use this definition in 
order to get a conservative upper bound on the BAU produced during the transition. 




(2) 
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For thick bubble walls the relevant baryogenesis mechanism is known as "spontaneous 
baryogenesis" , reviewed below. This mechanism, first introduced for baryogenesis during 



a second order phase transition [^5|, involves a space-time dependent field, which evolves 
coherently during the transition. This time evolution produces a CPT-violating term in 
the effective Hamiltonian called a "charge potential", which resembles a chemical poten- 
tial!. A charge potential will cause the free energy density inside the bubble walls to be 
minimized for nonzero baryon number, and hence the production of a net baryon num- 
ber via anomalous weak interactions. Spontaneous electroweak baryogenesis during the 



weak phase transition has been suggested for the two Higgs model in refs. |17,18| and as a 



baryogenesis mechanism for supersymmetric models in refs. [p0| , p!8[ . 

In the MSSM the ino mass matrices are space-time varying during the transition, and 
contain an irremovable CP violating phase, which leads to a charge potential for baryon 
number in the effective fermion Hamiltonian. If the inos are not too heavy and the phase is 
not too small, this charge potential is large enough to result in generation of an acceptable 
baryon number during the weak transition. We will find that the resultant BAU depends 
sensitively on the CP violating phase (pB, and on the masses of the inos . We are able 
to use the BAU to place upper bounds on ino masses and lower bounds on (f)B and the 
electric dipole moment of the neutron. 

Before launching into the specifics of the MSSM calculation, we show how to calculate 
the charge potentials resulting from space-time varying fermion mass matrices. 

Consider a fermion mass term of the form 

i)JCmij{x^)i)j + h.c. , (3) 

where we take all fermions to be left-handed, and C is the charge conjugation matrix . We 
can make a space-time dependent unitary change of basis on the fermions: 

tpi Uij{x^)ipj . (4) 



^ Recall that a chemical potential is a Lagrange multiplier introduced to implement a con- 
straint, which appears in the effective Hamiltonian as an energy splitting between particles and 
anti-particles. Similarly, a charge potential results in different energy levels for particles and anti- 
particles, but the energy difference is a real physical effect resulting from dynamical violation of 
CPT during a phase transition. 
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in order to make the fermion mass terms everywhere real, positive and diagonal; however 
the space-time dependence of U requires that we replace the kinetic energy terms in the 
Lagrangian by 

%.E. ^ %.E. + i^riuU^^U)^|; . (5) 
Note that since [/ is a unitary matrix, U'^ d^U may be written 

t/ta^t/ = za^5^a„(a;^)T„ , (6) 

a 

where for n fermions the T^s are generators of U{n), and the functions a{x^) are defined 
by eq. (|^). Thus the Lagrangian with mass term is equivalent to a Lagrangian with a 
real diagonal mass term but also containing a term 

- J]a^a„(x^)^T„7'^^ . (7) 

a 

If the transformation has a gauge anomaly there will also be a modification of the 
Lagrangian 

P I3a \ J 

where the tps are gauge generators in the left-handed fermion representation and the Fps 
are the gauge field strengths. 

The presence of these anomalous terms complicates the discussion of the charge po- 
tentials, and consequently we will choose our unitary transformation to have no gauge 
anomaly; this is the strategy followed in . 



Finally, if the transformation does not correspond to a symmetry of the interactions 
the coupling constants will be affected by the change of basis; however the effect of these 
coupling constants on energy levels is higher order in perturbation theory and will not 
concern us here. 

What effect does a charge potential have on a thermal system? There are two possi- 
bilities: 

1) If there is a charge potential for an exactly conserved charge, (e.g. electric charge or 
B-L) we can ignore it. Although it looks like the system could lower its free energy by 
producing a net charge, charge conservation imposes a zero charge constraint. It can 
easily be seen by integration by parts that a charge potential for an exactly conserved 



6 



charge has no physical effect; equivalently, we can always redefine fields so such charge 
potentials never arise. 

2) Charge potentials for non-conserved charges will lead to an asymmetry in the rates be- 
tween processes which create and destroy the charges, until the system reaches thermal 
equilibrium. The thermal occupation numbers in general will be different for particles 
and their CP conjugates, due to CP violation and the dynamical CPT violation from 
the space-time varying scalar fields. For a charge which is approximately (but not ex- 
actly) conserved the system will take a long time to reach equilibrium. If the system 
is near thermal equilibrium the net rate of charge production can be computed using 
thermodynamic arguments. For instance if there is a small charge potential as for the 
baryon number current and no other charge potentials!, for a system starting with no 
net quantum numbers the constraint of zero net baryon number can be implemented 
by introducing a baryon chemical potential 

fXB = -OlB ■ (9) 

The chemical potential /U^ is just the force of constraint on the system, i.e. the 
derivative of the free energy with respect to baryon number. Now the difference 
between the rates of anomalous processes which create and destroy baryon number is 
just proportional to the difference in the change in the free energy per event, leading 
to an anomalous baryon creation rate of 

PB = — = — 7^ — , (10) 

where F^ is the rate of anomalous baryon violating events per unit volume. The factor 
of 9 comes about because each anomalous event changes the free energy by 3fj,B and 
changes the baryon number by 3 units. 

We can now calculate the BAU produced in the MSSM during the weak phase tran- 
sition. We first compute the charged and neutral ino mass matrices, as functions of the 
scalar field expectation values, assuming equ. (|^), and find the transformation on the ino 
fields Ui{H) which renders their masses real, positive, and diagonal. Uj will in general be 
anomalous and would, by itself, give rise to an effective WW operator, as well as giving rise 



Note that anomalous weak baryon number violating processes can be affected by a potential 



for any charge generator whose trace over left handed fermion weak doublets is nonzero 
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to charge potentials for the inos via equ. (|^). By further making a space-time dependent 
baryon number rotation of the hght fermions we can remove this anomalous operator, at 
the cost of introducing a charge potential for baryon numbeiH: 

la.(/"'"^i^ivf/tm!^4.).S. (11) 

We then need to know how anomalous processes inside the bubble walls are affected by 
CP violating terms. Since the inos can quickly come into equilibrium with the ino charge 
potentials via ordinary non- anomalous interactions the ino charge potentials will not affect 
anomalous processes!. The main effect on anomalous baryogenesis will come from the term 
(|TTD. If the system is near thermal equilibrium then eqs. (|T0|) , and (^TJ) can be used to 
find the total baryon number density produced during the transitionB 



Pb = dtp 



We still need to know the rate of anomalous baryon density production Vb, which has 
been estimated in the symmetric phase to be 



Tb ~ «tkT* , (13) 

while in the broken phase it is vanishingly small. Unfortunately, there is currently no 
reliable way of computing Tb inside the wall where the scalar expectation values are 
changing. Furthermore it has been claimed, based on some 1 + 1 dimensional simulations 



This could also be seen simply by using the anomaly equation to replace FF with the 



divergence of the baryon current, and integrating by parts [20|. 

^ In the limit of vanishing gaugino mass or p, there is an additional approximate symmetry 
and the anomalous baryon production rate will be further suppressed. We always assume we are 
far from this limit, i.e. that the ino masses are not much smaller than the temperature. 

^ The effect of the spatial component of the charge potential will be to also produce a baryon 
number current inside the walls. We find this current has no effect on the baryon density produced 
in the thermal frame in the limit where the walls are thick. 
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||27|| that anomalous baryon production inside the waU is an inherently non-equilibrium 
process. (Dine suggests this is unlikely to be the case for a transition with thick walls 
and small latent heat such as occurs in the MSSM |]2^.) We will use a conservative 
estimate for the maximum baryon number produced by simply computing the integral 
using equ. (|T^ ) for r^. We think it is more likely that the baryon production inside 
the wall will be suppressed when the Higgs vevs become large. For instance McLerran 
29[] estimates that baryon production is shut off when the Higgs vev reaches about (1/2) 



its value inside the wall, while Dine et al. [ p^pO[ ] estimate that this shutoff occurs for vevs 
of order g^\JP/{ATi). In Figure 1 we plot the baryon number produced as a function of 
the value of the Higgs vev where this shutoff occurs for a typical choice of parameters; 
note that the dependence is approximately quadratic. We conclude that it is most likely 
that the actual baryon number produced during the transition will be between ~ 1/4 and 
~ 10~^ times smaller than our most favorable estimate. 

In order to find the allowed range of parameters in the MSSM, for each value of the 
gaugino and higgsino masses we adjust the CP violating phase 4>b to be as large as is 



compatible with current constraints on the EDMN [^. We then compute the BAU, from 
equs. (^2]) and (p!3|). The allowed range of ino masses are those for which the upper bound 
on the baryon to entropy ratio is greater than 0.4 x 10"^*^. For masses in the allowed 
range, we compute the minimum value oi (f)B which is consistent with the BAU, and use 



the calculation of Kizukuri and Oshimo in ref. to find the lower bound on the EDMN. 

In Figures 2a and 2b we plot the lower bound on the EDMN, as a function of and 
the wino mass parameter m2. The first figure has a phase cj)B near 0, while the second 
figure has the phase near tt. The bino mass parameter mi is assumed to satisfy the GUT 
relation mi = m2(5/3) tan^ 6w- The signs in fig. 2b indicate the sign of the EDMN, while 
in fig. 2a the EDMN is negativdl. The black excluded region corresponds to an EDMN 
greater than the current bound of 10~^^ e-cm pil], the dark grey region is an EDMN 



greater than 10~^^ e-cm, and the light grey region greater than 10~^^ e-cm. Thus for any 
portion of the allowed parameter space the EDMN is greater than 10~^^ e-cm. (Similarly, 
if we assume the selectron mass is equal to the squark mass, the electric dipole moment of 
the electron is greater than 0(3) x 10~^^ e-cm.) Throughout all of the region where the 
EDMN could be less than 10~^^ e-cm the lightest chargino mass is lighter than 88 GeV 



In computing the EDMN, we use the quark model and neglect the contribution of the phase 
')A', we have assumed the A-terms are small. 
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and the lightest neutrahno mass is hghter than 44 GeV. Also, if the EDMN is less than 
10~^^ e-cm the phase (pB must be near tt in order to avoid a chargino mass lighter than 
Mz/2. 

These figures have assumed that the rate of baryon violation is given by equ. |j2^, 



and that this violation is occurring throughout the bubble wall. If this is not the case, the 
baryon number produced will be reduced and the resulting phase must then be larger to 
compensate, giving rise to larger dipole moments. 

If the experimental bound on the EDMN is pushed down by an order of magnitude, 
or if numerical calculation finds the rate of anomalous baryon creation inside the phase 
boundary to be much less than it is in the symmetric phase, then there are strict upper 
bounds on ino masses. If the MSSM is responsible for baryogenesis the prospects for 
discovering electric dipole moments and supersymmetric particles in the next few years are 
excellent. Should we fail to make these discoveries, there are several possible conclusions. 
The BAU may come from physics above the weak scale, from weak scale physics other 
than supersymmetry, or from a more complicated supersymmetric model. The simplest 
extension of the MSSM would be to add a gauge singlet superfield, which changes the 
allowed form of the Higgs potential and greatly relaxes all the constraints on scalar masses 
|3'4 j. The tree level potential will generically contain large cubic terms, and gives a very 
strongly first order transition with a thin boundary between the two phases. Furthermore, 
with gauge singlets there is the possibility for CP violating phases in the Higgs potential. 
With thin bubble walls we expect baryogenesis to be dominated by the mechanism of 
refs. [1T6|JT9|] , in which CP violating particle scattering processes from the phase boundary 



leads to a transport of particle quantum numbers, biasing anomalous baryon production 
throughout the symmetric phase. The charge transport mechanism could produce the BAU 
for phases as small as 10~* ||19|| . Unfortunately such extended supersymmetric models 



currently have too many free parameters to allow for calculation of the BAU. 

In summary, we have shown that unlike the minimal standard model, the MSSM is 
still viable, but soon either it should be ruled out (for baryogenesis), or new experimental 
discoveries such as electric dipole moments will give us an important clue towards the 
explanation of why there are more baryons than anti-baryons. We are still a long way 
from testable predictions from baryogenesis in extensions of the MSSM. 
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Figure Captions 

Fig. 1. Baryon number density as a function of the Higgs vev. 

Fig. 2a. Lower bound on the EDMN for phase (f)B near 0. The black region is excluded, 
the grey region has an EDMN greater than 10~^^ e-cm, and the white region has 
an EDMN greater than 10"^'^ e-cm. 

Fig. 2b. Same as Figure 2a, except the phase 0b is near tt. The ± signs indicate the sign 
of the contribution of 0b to the down quark electric dipole moment at one loop, 
when 0B is chosen to produce a positive BAU. 
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